Singlet Superconductivity in Carbon Nanotubes 
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We discuss phase diagrams obtained by calculating temperature dependent correlation functions 
for superconductivity and density waves for a single (5,0) carbon nanotube (CNs). We use one- 
loop renormalization group method within logarithmic accuracy. In this system we must specify 
scattering channels in terms of momentum along the circumferential direction as well as the axis 
direction, since (5, 0) CN has two degenerate bands crossing the Fermi energy with circumferen- 
tial momenta. We find singlet-superconducting phase with 5?i circumferential angular momentum. 
Critical exponents for each correlation functions are also calculated. 

PACS numbers: 74.20.Mn 73.63.Fg 61.46.-(-w 05.10.Cc 



I. INTRODUCTION 

The discovery of carbon nanotubes (CN) in 1991 has 
attracted much attention because of its potential for new 
physics as well as applications in electronic devices • A 
single-wall CN (SWCN) is made of a graphite layer rolled 
up into a cylinder with small diameter. SWCNs may be 
regarded as one-dimensional (ID) electron systems due 
to quantization of the circumferential momentum, as con- 
firmed by experiments Q . 

Recently, Tang et al. reported evidence for super- 
conductivity in a SWCN system below 15K. This system 
consists of (5,0) CN with diameter 4A, separated from 
each other by zeolite walls. So this experiment implies 
that well-developed superconducting correlation (SCC) 
can exist in each individual nanotube. The microscopic 
origin of the SCC in a SWCN is an interesting ques- 
tion to ask, since divergent ID-SCC will eventually cause 
three-dimensional superconductivity with electron tun- 
neling through zeolite walls. 

On the other hand, bosonization |^ and renormal- 
ization group 0, 111 methods have been applied to the 
(n, n)CN to explain the superconductivity in the ropes 
of (n,n)CNs These works have shown that a sin- 
glet superconducting phase is more favored if short-range 
attractive interactions exist between electrons. More- 
over, some works ,ic| 0, have suggested the pres- 
ence of short-range attractive interactions in (n,n)CN. 
However, results of bosonization and renormalization 
group [ElESIll cannot be directly applied to (5, 0)CN 
case, because the band structure of (5, 0)CN is very dif- 
ferent from that of (n, ?t,)CN. 

Within the local-density approximationj_[5,0) CN has 
three bands crossing the Fermi energy |l,l |: one non- 
degenerate band and two two-fold-degenerate bands with 
opposite circumferential momenta, as shown in Fig.l. 



Such a band structure cannot be obtained within the 
tight-binding approximation because ct-tt hybridization 
effects change it significantly in a SWCN with diameter 
4A This band structure will provide opportunities 
to develop SCC with total circumferential momentum. 
We should also take into account the interaction between 
electrons with momentum transfer in the circumferential 
direction as well as in the axis direction. 

In this paper, we perform one-loop RG calculation and 
derive the scaling equations for the three-band system, 
including the interactions with angular momentum trans- 
fer. We find phase diagrams in terms of several coupling 
constants by solving these equations numerically. Criti- 
cal exponents which govern the temperature dependence 
of correlation functions near are also calculated. 



II. THREE-BAND MODEL AND THREE 
GROUPS 



We consider a system consisting of three linear bands 
near the Fermi level (Fig.l). Band has no angular mo- 
mentum (Lq = 0), band 1 and 2 have a finite angular 
momentum and are two-fold degenerate (Li^2 = i^'^i,2)- 
In general two-particle correlation (TPC) are formed be- 
tween branches with different sign in Fermi velocity in ID 
electron systems. As we have three bands, six channels 
of TPC are produced. We take forward and backward 
scatterings into account in this three-band system and 
perform g-ology 15], and determine the most divergent 
TPC among all possible ones by solving scaling equa- 
tions. We start from a Hamiltonian for non- interacting 
electrons: 
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where vf,'y and kp,j represent Fermi velocity and Fermi 
momentum for band 7. Spin, sign in Fermi momentum 
and in angular momentum L-y are denoted by cr, r and I = 
± respectively. We carry out a perturbation calculation 
with the interaction between electrons, 



Tl,T2,T3,r4 

rr \ ^ h,h,h,h 

^ / , ^71. 72, 73. 74 



7i.ri,h.-(^ + 9)^72,r2,i2,-'(^'-9) 

(k). (2) 



Keeping the momentum transfer in mind both in 
the circumferential direction as well as in the axial 
direction, we abbreviate interaction constants as gf 
where i and j represent the types of scattering pro- 
cesses in axial and circumferential direction respectively: 
1 =backward, 2 =interbranch forward, 3 =Umklapp and 
4 =intrabranch forward. We neglect the intra-branch 
forward scatterings f/4 because they are irrelevant within 
the one-loop RG framework. We also neglect Umklapp 
processes 53 because (5,0) CN systems are far from half- 
filling for any band. 

We divide six channels into three groups by considering 
the types of scattering processes. We denote a channel 
from band 7 and band 7' by (7, 7'). 

Group 1 : (0,0), (0,1), (0,2) 

These channels contain only two types of interactions 
(&o, 71/0,7) because of the restriction in momentum con- 
servation. Backscattering 69,7 exchanges the Fermi 
points of the incident electrons and forward scattering 
fo^j does not in contrast. 
Group 2 : (1,1), (2,2) 

These channels contain backward processes 52) the 
circumferential direction, which exchange angular mo- 
mentum between electrons. We will consider six inter- 
actions {gl,gh9i,9h9h92) here. 
Group 3 : (1,2) 

If \ni\ -f \n2\ = 5, these channels contain Umklapp pro- 
cesses {gi, ), which conserve the total momentum along 
the axial direction but does not conserve the total angu- 
lar momentum along the circumferential direction. We 
have to allow the angular momentum change by a re- 
ciprocal lattice vector. We will consider six interactions 
represented by {gl, gf, gf: g^: g^ 9t)- H \ni\ + \n2\ ^ 5, 
then we shall put gf = g2 = 0. 



III. SCALING EQUATIONS 

Our RG calculations are straightforward generalization 
of the method of Solyom |0| and performed by sum- 
ming up all one-loop diagrams i.e. Cooper-channels and 
Peierls-channels. 



A. Group 1 

These channels have two types of scattering processes 
(&; backward) and (/; forward). RG analysis of these 
channels is the same as pure ID system. RG equations 
for coupling constants are 



(60,7)' = 26^,^, 

(/0,7)' 



0,7' 



(3) 



Here 60,7 = 6o,7/27r?;o,7 and /o,7 — fo^-y /2t:vq^^ are nor- 
malized dimensionless coupling constants where wo,7 = 
+ vf,j)/2. The prime denotes differentiation with 
respect to x, where x = In with Ec being an energy 
cutoff. In these channels four types of TPC are produced. 
They are 



k,a=± 



d*-,(gr) ^ V'7,-,i,^(fc)V'o,+,.(gr - fc), (4) 



fe,CT=± 



where = kpfi i kp.-y represents the total momen- 
tum of the TPC. cdw, sdw, ssc and tsc stands for charge- 
density wave, spin-density wave, singlet supercontuctor 
and triplet superconductor, respectively. We can also de- 
rive scaling equations for these two-particle correlation 
functions x„i where Xm = (Ol^Om) and x„ 
and m stands for a kind of TPC. 





FIG. 1: Band structure of (5,0) CN. Electrons are labeled by 
band 7 £ {0,1,2}, Fermi point r = ±, angular momentum 
= ±h\nj\ and spin a — ±. 



dx 



logXm = K„ 



K 

K' 



0,7 

cdw 
0,7 



sdw 
^0,7 



0,7 

tsc 



4&0,7 — 2/0,7, 

^2/0,7, 

26o,7 + 2/0,7, 
— 26o,7 + 2/0,7. 



(5) 



We determine the most divergent TPC by solving these 
scaling equations. Phase diagram for group 1 is the same 
as that of a pure ID system with backward and forward 
scatterings as shown in Fig. 2. When we have 60,7 < 
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FIG. 2: Phase diagram for group 1 in (/o.7- plane. Phase 
boundary is given by two lines 60,7 = 2/0,7 bo,j = 0. 

initially, it is scaled towards negative large quantity 
and diverges at finite temperature. We call this region 
gapful because gap opens in tlie excitation spectrum at 
tliis temperature. 

B. Group 2 

These channels have six scattering processes shown in 
Fig. 3. Scahng equations for these six couplings are 



i~9iy = 




+ 2glgl, 




(4.9f + 


2~gl - 2~gt)~gl + {2gl - 2gf)g 


(gtr = 




+ 2{gt)^-2gf~gl, 


{-gir = 




+ 2{gl-~g^)~g\, 


{9I)' = 


igir - 


^{gff + glgl, 


{gty = 


i-gfr - 





where gf = gj 12-kvf,-^ are normalized coupling constants. 
Couplings are scaled towards large values and diverge at 
some finite temperature except for the case of g\—g\ = 
0. This suggests that the channel from group 2 tends to 
produce gapful phase when it contains interactions with 
angular momentum transfer. 

There arc several types of TPC in these channels and 
these orders have either zero or ±2n-y total angular mo- 
mentum. Symmetry in spin space and in circumferential 
direction characterize each TPC; either of symmetric or 
antisymmetric order with respect to inversion will appear 
both in the spin and the circumferential direction. We 
denote all the twelve types of TPC as; 



k,cr,a' 

Hi, 1^2 = ±, L7 = fin^, o-±, ^± = ^ Q ± 1 ^ ' ^±^M±= Q ^ . 
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We have abbreviated density wave and superconductor 
as dw and sc in the superscript, fxi represents symmetry 
in spin space and fj,2 represents symmetry with respect 
to circumferential direction. We use + for symmetric 
orders and — for antisymmetric orders. Therefore fii — + 
means CDW and TSC and m = — represents SDW and 
SSC. Scaling equations for these two-particle correlation 
functions have the same form as the first line in equations 
(5). Ki for each orders can be calculated as 

Ki^(O) = (4.9^ - 2.92*) + - 2gl), 

K''''"i2ln^) = Agl - 2gl 
K^'''^{2ln^) = -2gl 



K^^%2ln^)^2~gt + 2~gl 

K''%2ln^) = -2~gl + 2~gl 

^r.'^(O) = (25i + 2g) + M2(25f + 2~gl), 

^^T(0) = (-251' + 25^) + M2(-25? + 2gl). (8) 

By substituting the solution for equations (6) into (8), 
we obtain the temperature dependences of two-particle 
correlation functions. We then obtain the phase diagram 
in coupling space by specifying the most divergent corre- 
lation function. The phase diagrams are shown in Fig. 4. 
If we assume bare value for couplings which depends only 
on the momentum transfer, we can put "5^ = gf = b, 
ql = 92 = /" and get four independent parameters 
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For the special case "gj = 92 = 0", we can solve the 
equation (6) exactly, and the solution is 



&(0) 



b{x) = 

l-46(0)x 

b(x) - 4/(x) = 6(0) - 4/(0). 



(9) 



When &(0) < 0, & diverges at finite temperature Tc = 
i?c exp(l/4&(0)). In one dimension, such an anomaly at 
finite temperature is an artifact of one-loop calculation, 
and higher-order terms will shift it to T = 0. Never- 
theless, this Tc suggests us the temperature where the 
system crosses over from weak coupling to strong cou- 
pling. 

The phase diagram for this special case is shown in 
Fig. 4(a). In gapless phase, all the DWs and SCs are 
degenerate and have the same temperature dependence 
of the correlation function. The DW phase and the SC 
phase are separated by a phase boundary 5 — 4/ = 0. In 
the gapful phase, only CDW+(0) phase appears, where 
+ and (0) represent the symmetry /i2 = + and angular 
momentum L = respectively. 

In general, gl and take finite values, and other 
gapful phases enter into the phase diagram as shown in 
Fig. 4(b). The effect of / on the phase diagram is negligi- 
ble. Independent of the values of b and /, there exist four 
gapful phases in which the most divergent order of each 
is SSC+(0), SSC_(0), CDW+(0) and CDW_(0) respec- 
tively. There is no gapless phase except for the original 
point (gljCjl) — (0,0) even when b > 0. Four gapful 
phases are characterized by asymptotic behavior of the 
couplings near the critical temperature Tc [Ifij . which 
are given by 



/ = + 
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FIG. 3: Six scattering processes in channels of group 2 
are shown in momentum space when (no,n\,n2) — (0,3,2). 
Hexagon represents the Brillouin zone. 



gi = gtA/il + Ax] 
Tc = Ecexp{-A-^) 



(10) 



All the couplings are proportional to A/[l + Aa:], which 
diverges at Tc given above. Sets of gj* are universal for 
each phase as in , and listed in Table ^ Equations 
(6) are invariant to transformation (ffi , 52)~^(~ffi j "ffD- 
The two sets of solutions (gl* , gf* , gf* , g2* , 9^* , gt*) and 
{9l*,-9l*,gt*,-9l*,g^*,92*) correspond to two phases 
with different symmetry {^2 = ±) as seen in Table 
The asymptotic solutions lead to behavior of TPC, 



Xni oc I In 



|a™-l-l 



Tc 



oc 



'T-T 



(11) 



where am determines the critical exponent for the corre- 
lation function of TPC of type m and is calculated from 
the value of gj* using a„i — Kmlgj^gj' ■ Correlation func- 
tions diverge at Tc if the exponents are positive, and are 
suppressed if they are negative. These exponents in each 
phase are listed in Table HTl 

It can be seen from Table |2 that scattering processes 
with angular momentum transfer gl and are renormal- 
ized to in the two CDW phases, on the other hand they 
are renormalized to large value in the two SSC phases. 
This means such scattering processes become strong at- 
traction between electrons in a SC pair and give rise to 
divergence of SCC. We should note that all the dominant 
TPC in these channels carry zero angular momentum, 
though they have the same correlation function as TPC 
with angular momentum ±2hn^ in the non-interacting 
case. Only one type of TPC diverges and all others are 
suppressed in every phase. 

We see from numerical calculations that Tc expo- 
nentially depends on backscatterings b and Tc becomes 
higher with negative large b. However Tc will quickly 
be suppressed in case of positive b. Therefore phase 
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FIG. 4: (a) Phase diagram for group 2 in (/, b) plane is 
exactly obtained for the special case gl = gl ~ 0. (b) Phase 
diagram for group 2 in (^1,52) plane. Fixed (6,/) values 
are shown above each diagram. In the shaded region |T^ . 
we cannot determine which TPC is the most divergent since 
couplings diverge at very low temperature. 
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transition can occur at 15K in these channels only when 
the backscatterings b are attracting forces. However, T^. 
tends to be higher than that of group I's for two rea- 
sons; (1) double degeneracy of opposite angular momenta 
and (2) smaller curvature of the bands ~ vf,2 = 
2.8 X 10^ m/s < VF.a = 6.9 x 10^ m/s) will produce a 
larger density of states of electrons compared to that for 
band 1. 



C. Group 3 

This channel has six independent scattering processes 
shown in Fig. 5. Scaling equations for these six couplings 
are 

{-glY = 2{~gir-2~gfgl 

{gfy = {4gl+2~g^^^2gM + i2gt-2gl)gl 
{-gty = 2{gtr + 2~g!gl 

{glY = {9\f-{9l)\ 

{gl)' = 2gfgf + {2gl-2gl)gl 

i9iy = {9ir+{9tr+i9ir, (12) 

where gf = gl /'n'{vF,i+VF,2) in (10) are normalized scat- 
tering amplitudes. All couplings are scaled to large val- 
ues except for the case = = 0- This suggests that 
the channel from group 3 tends to produce gapful phase 



when it contains Umklapp processes g\.,g2- 



l = + 




FIG. 5: Six scattering processes in channels of group 3 
arc shown in momentum space when (no,ni,n2) = (0,3,2). 
Hexagon represents the Brillouin zone. 

There are twelve types of TPC: 



k,cr,a',l,l' 

or.(93-,L+) = 4^ E V'2,-,i,.(93- - k) {dlf ® M^f ) Vi,+,;',.'(fc), 

OZ,^M^l^~) = 4f E V'2,-,-;,.(?3" - k) dlf i^i,+M'ik)- (13) 



We use symbols cr^i, (t^i, M^2 and Mfj_2 which arc the 
same as the ones in (7). TPC carry momentum in the 
axis direction qf = kp,! ^ kp,!- They carry angular mo- 
mentum L"*" or L~, where L"*" = h{ni + ^2) = 5h and 
= h(ni — 712). The Umklapp processes make transi- 
tion between TPC with -L+ and TPC with L+ in the 
first order correction, therefore bonding or anti-bonding 
states are formed. In other words, gf and §2 produce the 
symmetric or anti-symmetric states with respect to the 
inversion of circumferential direction. 

Ki in the scaling equations for these two-particle cor- 



relation functions are given by 

(L+) = (45i' - 252') + M2(45? - 2M), 

<"(L+) = -2(5i + M25i), 

if-;(L+) = {2gt + 2gl) + ^2(2^? + 2gl), 

= i~2~gt + 2~gt) + M2(-2ff? + 2gll 

K^^\h-)^2g\ + 2~gl 

if*-(L-) = -2ffi + 2ff2. (14) 
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TABLE I: Strong-coupling fixed-point parameters for eacli gapful phase in group 2 



phase 


gl* 


9? 


gt 




g^ 


gl' 


CDW+(0) 





-1/4 


-1/4 





-1/16 


-1/16 


CDW-(O) 





1/4 


-1/4 





-1/16 


-1/16 


SSC+(0) 


-1/4 


-1/4 





-1/4 


-3/16 


1/16 


SSC-(O) 


-1/4 


1/4 





1/4 


-3/16 


1/16 



TABLE II: Exponents -(a™ -I- 1) in group 2 



phase 


CDW+(0) 


CDW-(O) 


CDW(2Lt,) 


SDW+(0) 


SDW_(0) 


SDW(2Lt,) 




SSC+(0) 


SSC-(O) 


SSC(2Lt,) 


TSC+(0) 


TSC-(O) 


TSC(2Lt,) 


CDW+(0) 


7/8 


-9/8 


-9/8 


-9/8 


-9/8 


-9/8 




-3/8 


-11/8 


-3/8 


-11/8 


-3/8 


-11/8 


CDW-(O) 


-9/8 


7/8 


-9/8 


-9/8 


-9/8 


-9/8 




-11/8 


-3/8 


-3/8 


-3/8 


-11/8 


-11/8 


SSC+(0) 


-3/8 


-11/8 


-3/8 


-11/8 


-3/8 


-11/8 




7/8 


-9/8 


-9/8 


-9/8 


-9/8 


-9/8 


SSC-(O) 


-11/8 


-3/8 


-3/8 


-3/8 


-11/8 


-11/8 




-9/8 


7/8 


-9/8 


-9/8 


-9/8 


-9/8 



If we have gf = = initially, these equations can be 
solved and phase diagram becomes equivalent to that of 
group 1. In this case scaling equations (12) are reduced 
to 



(.9i)' = 



2(5^)^ 



gl and gf diverge at different temperature Tc = 
i?c exp(l/2g}) and = i?c exp(l/2,gf) respectively. 
Gapless phase appears in the region gl,gf > 0, and gap- 
ful phase appears in the region gl,gf < 0. 

By solving (12) and (14) numerically, we obtain phase 
diagrams (Fig. 6) for the case where Umklapp processes 
have some finite value (gi,52 0)- 



(-0.06,-0.02) 



CDW,(L) 


SSC_(U) 






SSC,(L-) 


cdw-(l:) 


0.2 -0 1 ( 


0.1 _ 


(0.02,-0.02) 


cDw.(i;) 


SSC.{IJ) 


SSC,(L') 


COW-(L) 



(-0.02,-0.02) 



CDW,(U) 



SSC_(L) 



SSC,(L') 



CDW-{L) 



.0 2 ^1 Q. I 



(0.06,-0.02) 



-0.2 -0.1 a 0.1 2 



CDW.(C) 


SSC.{LJ ) 




SSCi(L') 


CDW-CL) 



We have taken gl = gf = b and = .92 ~ / f^'' ^i^^^" 
plicity. The effect of / on the phase diagram is negligible. 
Independent of 6, / values, phase diagram mainly consists 
of four large gapful region {CDW±(L+), SSC±(L+)}, 
and there is no gapless phase when we have non-zero 
Umklapp processes (gi,^!)- These gapful phases are 
characterized by the asymptotic solutions (10) and (11) 
as well as in group 2. We show the values of gf* in Table 
mil and exponents — (ofm + l) in Table llVl Equations (12) 
are invariant to transformation (ffi , 32)~*(~5i j "Sl)- 
The two solutions (ffi*, fff , fff, fff, sl*, sf) and 
(5'i*,-5i*,5f,5l*,-52*,52*) correspond to two phases 
with different symmetry (/i2 = ±) as seen from Table 

EH 

In CDW^2 (L^) phase, the correlation function for both 
CDW^2 (L"*") and TSC-^^ (L+) diverges at and all oth- 
ers are suppressed. Correlation function for SSC^2(L"'") 
alone diverges in SSC^2(L~'') phase. Every divergent 
TPC carry angular momentum = 5h. Umklapp pro- 
cesses (^1 , gl) are scaled towards large value in any case, 
which implies that they play an important role for TPC 
with angular momentum to develop. Tc is exponen- 
tially dependent on the value of backscatterings b as in 
the case of group 2. When 6 > Tc is suppressed to very 
small energy scale. This seems to be a trace of gapless 
phase in the special case gf — g2 — 0- When b < 0, 
however, Tc will be higher than that of the gapful phase 
in group 1 for the same reason stated in group 2. 



-0,2 1 0. 1 2 



FIG. 6: Phase diagram for group 3 in (ffijfff) plane. Fixed 
{b, /) values are shown above each diagram. In the shaded 
region, we cannot determine which TPC is the most divergent 
since couplings diverge only at very low temperatures, where 
our calculation becomes less reliable. 



IV. SUMMARY 

We have discussed possible phase transitions in (5,0) 
CNs, which is quasi-lD system with angular momentum, 
by evaluating the most divergent two-particle correlation 
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TABLE III: Strong coupling fixed point parameters for eacli gapful pliase in group 3 



phase 


or 


gr 


gt 


gf 


gr 


gt 


CDW+(L+) 


-0.1258 


-0.4219 


0.0810 


-0.0034 


0.1116 


-0.1970 


CDW-(L+) 


-0.1258 


0.4219 


-0.0810 


-0.0034 


-0.1116 


-0.1970 


SSC+(L+) 


0.1026 


-0.2302 


-0.2764 


0.0616 


-0.2685 


-0.2015 


SSC_(L+) 


0.1026 


0.2302 


-0.2764 


0.0616 


0.2685 


-0.2015 



TABLE IV: Exponents -(q^ + 1) in group 3 



phase 


CDW+(L+) 
SSC+(L+) 


CDW-(L+) 
SSC-(L+) 


CDW(L-) 
SSC(L-) 


SDW+(L+) 
TSC+(L+) 


SDW-(L+) 
TSC-(L+) 


SDW(L-) 
TSC(L-) 


CDW+(L+) 


1.4073 
-0.1474 


-2.4143 
-1.3886 


-1.7182 
-0.7416 


-0.7836 
-1.5109 


-1.2299 
0.6231 


-1.3940 
-1.2449 


CDW_(L+) 


-2.4143 
-1.3886 


1.4073 
-0.1474 


-1.7182 
-0.7416 


-1.2299 
0.6231 


-0.7836 
-1.5109 


-1.3940 
-1.2449 


SSC+(L+) 


-0.9034 
0.9531 


-1.6708 
-1.0416 


-0.2973 
-1.3283 


-1.4139 
-1.0732 


-0.3398 
-1.2265 


-1.4030 
-0.9180 


SSC_(L+) 


-1.6708 
-1.0416 


-0.9034 
0.9531 


-0.2973 
-1.3283 


-0.3398 
-1.2265 


-1.4139 
-1.0732 


-1.4030 
-0.9180 



function. 

We have found new type of superconducting phase, 
where a Cooper pair carries non-zero angular momen- 
tum and Umklapp processes play significant role. All 
other dominant TPC are singlet-superconducting or 
charge-density wave with zero or angular momentum, 
which are caused by backscattering or Umklapp scat- 
tering. Back scattering in the circumferential direction 
strengthen the correlation between electrons with oppo- 
site angular momentum and cause divergent TPC with 
zero angular momentum. On the other hand Umklapp 
processes strengthen the correlation between electrons 
with angular momentum in the same direction, and cause 
divergent TPC with angular momentum. 

Tc evaluated from RG equations become quite low with 
repulsive interactions between electrons. So we need a 
large reduction of Coulomb interactions or strong attrac- 
tions by electron-phonon interactions for the phase tran- 
sition to occur at 15K. However, in group 2 and 3, Tc will 
be higher than that of group 1 because they have higher 
density of states at Fermi level. Therefore it seems that 
group 2 and group 3 play a critical role in the supercon- 
ductivity observed in (5,0) CN. A recent work [iJI has 
shown that acoustic phonon exchanges result in strong 



attractive interactions in nanotubes with smaller diame- 
ter. Thus, there may be some possibility of understand- 
ing the unusually high transition temperature (15K) with 
our result. 

In our study we have not considered the effects of 
pair-tunneling processes such as {(7, ±r), (7, =Fr)} 
{(7', ±t), (7', Tt)}- Such processes are known to be im- 
portant for the formation of SCC in (n, n) CNs H 0, Q 
and in two- and three-leg ladder systems ^ well. 
Therefore it is also important to study the effects of pair- 
tunneling processes. However things will be more difficult 
because pair-tunneling processes will mix with the TPC 
in every channels and we cannot treat three groups in- 
dependently. Effects of pair-tunneling processes will be 
discussed elsewhere. 
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